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Logocyclic Curve by a continuous motion ; and a very ingenious 
instrument has been contrived by Mr. Henry Johnson of Crutched 
Friars, to describe the spiral of Archimedes, which is as simple as 
it is effective. 


June 17, 1858. 

The LORD WROTTESLEY, President, in the Chair. 

The Earl Granville, Professor Hennessy, and the Rev. Samuel 
Iiaughton were admitted into the Society. 

In accordance with notice given at the last Meeting, the Earl of 
Rosse proposed the Right Hon. Sir John Pakington, Rart. for 
election and immediate ballot. 

The Ballot having been taken, Sir John Pakington was declared 
duly elected. 

The following communications were read 

I. “ On the Problem of Three Bodies.” By Charles James 
Hargreave, LL.D., F.R.S, Received May 3, 1858. 

(Abstract.) 

The author states that the principal object of this memoir is to 
set forth two new methods of treating the dynamical equations by 
processes of variation of elements, differing from the ordinary pro¬ 
cesses of this nature principally in this particular, that the variations 
are represented in explicit terms of the elements themselves and of 
the time, and not through the medium of partial differential coeffi¬ 
cients. It has been his object to render the processes as elementary 
as possible; and to preserve them in a vigorous form, by post¬ 
poning all attempts at approximation until the formulae are actually 
applied to practical problems. The applications given, in the paper 
comprise the circular and spherical pendulums, and the planetary 
and lunar theories, and a special theorem as to the movement of 
the plane of a planet’s motion under the influence of several other 
planets. 

The original normal problem which is taken as the basis, is that 
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of motion about a fixed centre of force* where the force is directly as 
the distance; or, in other words, the system of equations not ex¬ 
ceeding three in number, of the form 

x ff 4 - n 2 x~0 ; 

whose solutions are represented under the form 

%=z\ a a cos (nt+p) b sin (nt + p ), 
y=\ b a cos (nt-\-p) -{-fi b b sin (nt+p), 
z=:\ e acos (nt+p) ~h/u e b sin (nt + p ); 

where 

cos (j) cos \p — sin <£ sin \f/ cos Z, 

\ b = cos (j> sin \f; -f sin (j> cos cos Z, 

A c = sin <p sin Z ; 

— sin <£ cos ^ — cos 0 sin ^ cos Z, 
yu^= — sin </) sin \p + cos </> cos \p cos Z, 
p c = cos (j) sin Z; 

to which are afterwards added, 

v a ~ sin d> sin Z, 

3 /^= —cos if/ sin Z, 

r c = cos Z. 

These are the equations of an ellipse whose centre is at the force, 
and situated in a plane inclined at the angle Z to the plane of x y , and 
the longitude of whose node is \p ; and 0 is the angular distance of 
the major axis of the ellipse from the node; a and b are the semi¬ 
axes of the ellipse; and p is the angular distance, from the major 
axis, of the zero-point of the motion, measured on the circle described 
on the major axis. A uniform motion around the circle represents 
the place of the body by the corresponding point on the ellipse, 
where it is cut by a perpendicular dropped on the major axis. 

If the force be not situated at the origin, but at the point (X, Y, Z), 
we have merely to substitute a?—X for x , &c. in the above equations 
of motion and solutions. 

Applying the method of tangential variation to the system 

x tf -\-ArX=0, &c., 
r s 

we perceive that this system admits of complete solution in finite 
terms, leading in fact to the usual theory of elliptical motion. Taking 
this system, therefore, as a normal system, the author proceeds to 
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deduce the formulae for the variation of the elements of this system, 
in order to arrive at the solution of the system 


+&c. 
r 3 


The elements which have been selected, for reasons fully explained 
in the paper, are l and \p, whose meanings are already known; A 
and Nr denoting respectively the mean distance, and the longitude 
of the epoch measured in the plane of the tangential ellipse as it 
exists at the time t, and measured from the node at that time; and 
e and w denoting respectively the eccentricity of the tangential 
ellipse, and the longitude of its perihelion measured as above ; and 
it is observed that these are strictly normal elements, according to 
Professor Donkin’s definition of normal elements. 

The variations of these elements are then rigorously found, and 
are expressed as follows :—Denote 

cos \p P A . + sin \p P^ by the symbol P^, 

and 


cos l (cos ^P*,—sin \pP & ) -f V z sin l by the symbol P,, ; 

and let 

— P^ sin 0 + P^ cos 6 == P^ 0 ; — P^ sin •ur -f P^ cos = P^ w ; 
Pf cos 0 -f P, sin 0 = P,^ e ; P^ cossr + P^ sin ra'—P-^; 

then 




r cos 0 


*(*.). 


~NA 2 (1 —e 2 )^ v 
r sin 0 

sm^^=NA X i_ e2)i (P0> 


ce= NA 2 (] { (c°s(0-sr) + tf)P fj9 +(l + ecos(0-w))P f> 

<?(J® + cos l Sf) = | sin (0 - w) F ( ,« -( 1 + e cos(0-w)) P„ >w | , 

2 

^—N (1 — e 2 )^(P^’'° e P^’ 

S(N(i+r))=^p^|—ffln(0—»)Pfc«+(l + ecos(0—w))P„ >M —2eP,. 


which are capable of being expanded in terms of the elements, and t 
by means of the ordinary expressions for r, 0 , and 0—w in terms of 
the same quantities. The values of the elements at the time t being 
supposed to be found, by the integration of these formulae, in terms 
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of t , and their initial values, are to be substituted in the ordinary- 
expressions for the coordinates, so as to obtain their values at the 
time t . 

The author exhibits the application of the preceding formulae to 
certain simple examples, and then proceeds to apply them to the 
planetary theory. For two planets (distinguished by the suffixes 2 
and 3) supposed to move in the same plane, the following are the 
rigorous expressions for the variations. Let a 2 and a 3 be the ratio 
of the mass of each planet to that of the central body. Let P denote 
the cube of r 3 -^r 23 , and let (P—1) sin (0 3 — d 2 ) be called Q, and 

(P—1) cos (d 3 —0 2 )—— P be called R; then 

r s 

{ ( cos ( 2 + ^ cos (0 2 —OHO 

Q + sin (0 2 —rav) (1 + e 2 cos (0 2 —0)R j , 

ejw a =et 3 { sin ( 2 + e * 

e 2/ L 

Q —cos (0 2 — & 2 )(l+e 2 cos (0 2 —nr 2 ))Rj , 
SA 2 =2a 3 a ^ aA % T i { (1 + cos(0 2 cr 2 )) Q+e 2 sin (0 2 -® 2 )R } , 

n( l ~ e D 1 > 

<5(N(£+t-)) = - (1 - e 2 )* hr - 2a 3 N 4^R. 

? 3 

From these formulae, the secular variations of the elements are 
obtained without difficulty ; and a new method of integrating the 
equations for the variations of the eccentricity and longitude of 
perihelion is given. 

The author then enters upon a minute examination of the mathe¬ 
matical character of secular variations, and their bearing Upon the 
methods of approximation to which the problem of three bodies 
has given rise. It is pointed out that the disturbance finally effected 
through the medium of a secular variation is not of the order of the 
disturbing force, or rather of the ratio of the disturbing force to the 
central force; but that it may remain precisely the same, though this 
ratio should be diminished or increased without limit. The differ¬ 
ence affects not the aggregate amount of deviation or disturbance 
caused, but the time in which this aggregate amount is produced. 

If we consider the undisturbed problem of two planets about a sun 
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as representing motion in two planes inclined to each other at the 
angle I, and in ellipses having eccentricities e 2 and e 3 , it is shown 
that, no matter how small or how large may be the disturbing force 
produced on each orbit by the other planet, the aggregate amount of 
disturbance of the planet m 2 is of the order of the quantities I and e 2) 
and that of the planet m 3 , of the order of I and e s . From considera¬ 
tions of this nature, which are dwelt upon at length in the memoir, 
the author concludes that the ordinary direct methods of solution 
by approximation, being based upon the erroneous assumption that 
the variations of the coordinates are of the order of the disturbing 
force, are not, in a mathematical sense, legitimate processes ; and 
that, in the planetary theory, they produce results practically true 
only on account of the minuteness of the disturbing forces, and the 
consequent great length of the secular periods ; and that, in the 
lunar theory, their failure is made evident, in consequence of the 
comparatively large magnitude of the disturbing force, and the con¬ 
sequent rapidity with which the elements of the moon’s orbit pass 
through their secular periods. 

The formulae for the variations of the elements are then applied 
to the lunar theory; and some of the integrations are effected 
by means of a lemma containing the solution of the differential 
equation 

d(f) . 

~^-=F cos (pt—q<f) 

(where F, p and q are numerical coefficients), in the form 


. x tfFcos Mi n 2 F 2 )^ 

cos (pt-q<p )= p+?F - oiM7 , M being (1- ^)p. 


By this method, the total motion of the moon’s perigee, as well as 
the coefficients of the evection, are fully obtained in the first in¬ 
stance, without the necessity of any second approximation; and the 
usual difficulty as to the movement of the perigee does not present 
itself. The motion of the node, and the evection in latitude, are 
correctly obtained in a similar manner. 

This part of the memoir is concluded by an extension of the 
general formulae for the tangential variation of elements to the case 
in which we suppose the constant fx to become variable, the result 
being to add to each variation a term involving fyi. 

The third part of the Paper contains the development of the 
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method of osculating variation, before briefly described ; from which 
are deduced the formulae for the osculating variations of elliptic ele¬ 
ments. This method is capable of being applied to the planetary 
and lunar theories, as well as that of tangential variation; but the 
advantages of this method did not appear to be such as to justify 
the actual expansion of the formulae for these theories. The author, 
however, shows that with reference to any system of three bodies, 
the equations of motion for each body naturally assume the form 

—X) = 0, &c. 

(being the system solved by this method) ; and that the X, Y, and 
Z are absolutely the same for each of the three bodies. This is 
shown by demonstrating, that at any given moment the three lines 
which represent the direction of the force acting on each of the 
three bodies all pass through the same point, which is denominated 
the centre of force. The coordinates of this common centre of force 
are, 

(23)*, + (31)* a +.(12)* a . 

(23) +(31) +12 ’ 

with similar expressions for Y and Z; (12) being r 12 ~mm 2 xr 129 
X denoting the law of force, &c. Each body has its own value of 
n 2 ; their ratios being denoted by the proportion 

n\ : ’n\ : n\ : : 

X r 23 X r n X r ™ 

The invariable plane of this system of three bodies is then found; 
and it is shown that the nodes of the three orbits upon this plane 
are always in a certain relative position, constituting a kind of tri¬ 
angle of equilibrium about the centre of force; resulting, in the 
limiting case where one of the three bodies is infinitely larger than 
the other two (or in what is denominated the undisturbed Problem 
of Three Bodies), in an exact opposition of the two nodes of the 
orbits of the latter two bodies upon the invariable plane of the 
system. 

The formulae for the osculating variation of elements are then 
applied to a system of three bodies, of which one possesses a pre¬ 
dominating magnitude, so far as is necessary to determine the move¬ 
ment of the planes of the orbits ; and it is readily shown that, if we 
consider only the first order of the disturbing force, the inclination 
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of the plane of each orbit to the invariable plane is absolutely con¬ 
stant ; and that the two nodes are always in opposition to each 
other, and move with a uniform angular velocity round the inva¬ 
riable plane. 

This theorem is then extended to a system of n bodies moving 
about a central predominant body; and it is shown that the ag¬ 
gregate effect of the disturbing forces of such a system upon the 
plane of any one of the bodies can always be represented by stating 
that its node upon a certain fixed plane revolves with a uniform 
angular velocity, the plane of the orbit always remaining at the 
same inclination to the fixed plane. The rate of this angular move¬ 
ment, and the coordinates of the fixed plane upon which the move¬ 
ment takes place, are found by means of formulae of remarkable 
simplicity. These three quantities may be ascertained once for all 
for each planet (viz. the inclination of the fixed plane on which the 
node moves to any coordinate plane, the longitude of the node of the 
fixed plane in relation to any coordinate line, and the angular rate of 
movement of the node of the orbit upon this fixed plane), and, when 
once ascertained, may be regarded as fixed elements of the planet, 
from which the position of the plane of its orbit can always be deter¬ 
mined without the use of tables. 

It is then shown that a system of the form 

x ,f -\-n 2 a ?=P#, &c., 

where n 2 and P*, P^, and P* are any variables, may be solved by the 
same set of final integrals, and the same values of x\ y\ and 
by supposing the elements a, b , </>, \p, z, l , and p to become variable. 
These elements are those of an ellipse tangential to the actual 
curve of motion ; and the following formulae are obtained for their 
variation 
Let 

f l a Pa’4- l l b P# 4“Pc P^ — (P#)> 

^^+>’4 Py + |/ cP^—(P*) ; 

and let (putting T for nt-\-p) 

a cos (j) cos T — b sin <f> sin T=£, 
a sin <p cos T + b cos 0 sin T= rj ; 

VOL. IX. V 
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then 

nab d l=£(P g ), 
nab sin l ^=??(P*), 

£(nab) = a cos T(P^) — b sin T(P^), 

\ln 2 (a 2 + b 2 ) — —n(a sin T(P A ,) — b cos T(P y )) + r 2 nn!, 

1 . ... rjd 

(a 2 ~ b 2 ) (^0+ cos l (hp)=—(b cos T(P#) — a sin T(P 2/ )) + 2a b sin T cosT—, 

1 n ! 

(a 2 —b 2 ) %(nt + p) = ——(a cos T(P A? ) — b sin T(P y))~(a 2 +b 2 ) sin T cos T~* 

It may be observed that £ and rj are coordinates of the body referred 
to the plane of the tangential ellipse, and to an axis of £ coinciding 
with the node. 

This method is denominated the method of Tangential Variation ; 
and it is applied directly to the problem of the circular pendulum, 
that of the spherical pendulum, and that of the motion of a particle 
where the force is a function of the distance, and in particular that 
of elliptical motion, where the law of force is that of the inverse 
square. 

In a subsequent part of the paper it is shown that a system of 
the form 

x 11 + ri 2 (x —X)=0, &c., 

where n 2 , X, Y, and Z are any variables, may be solved by the same 
set of final integrals, and the same values of y\ and z f as those 
which have been already given as the solutions of the same system 
when n, X, Y, and Z are constant, by supposing the elements to 
become variable. In such a case, the elements are those of an 
ellipse osculating with the actual curve of motion, always of course 
having its centre at the moveable point (X Y Z). The following 
formulae are obtained for the variation of these elements :— 

Let 

XaX' + X*Y' + XeZ' = (X'), 

p«X'H-^Y' + p,Z'=(Y'), 

VaX' + i'ftY'+ *£'=(#); 


fit 

nab 

nab sin l — 


then 
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S(nab)z= -»((X') b cosT + (Y') a sinT), 

P(n 2 (« 2 + b 2 ))= -k 2 ((X') a cos T + (Y') b sin T) + r'nri, 

(a 2 — b 2 ) (ty + cos l ty)= -((X') b sin T + (Y') a cos T) + 2 ab sin T cos T—, 


(a 2 -b 2 ) S(nt+p)=(X') asinT + (Y')6cosT-(a 2 + 6 2 )sinTcosT^-; 

in which ~ and — are the differential coefficients of the expressions 

for £ and r\, taken explicitly with regard to t. 

This method is denominated the method of Osculating Variation. 


II. “Description of some Remains of a Gigantic Land-Lizard 
{Megalania prisca, Ow.) from Australia.” By Prof. 
Richard Owen, P.R.S. Received May 13, 1858. 

(Abstract.) 

The subject of this communication forms part of a collection of 
fossil remains from Australia, recently acquired by the British Mu¬ 
seum, and demonstrates the former existence in that continent of a 
land-lizard considerably surpassing in bulk the largest species now 
known. The characters are chiefly derived from vertebrae, partially 
fossilized, equalling in size those of the largest existing Crocodiles ; 
they are of the c procoelian , type, hut present lacertian modifica¬ 
tions, and closely agree with those in the great existing c Lace-li¬ 
zard ’ of Australia (Ilydrosaurus giganteus , Gray), of which indi¬ 
viduals upwards of six feet long have been taken. A generic or sub¬ 
generic distinction is indicated by the comparatively contracted area 
of the neural canal, and by the inferior development of the neural 
spine, of the fossil vertebrae, which have belonged to an indivi¬ 
dual not less than twenty feet in length, calculated from the vertebrae 
and proportions of the body of the existing Hydrosauri . For this, 
probably extinct lizard, the name of Megalania prisca is proposed. 

The results of an extended series of comparisons of its vertebrae 
with those of recent and extinct Sauria are given; and the paper is 
illustrated by drawings of the vertebrae of Megalania and those of 
Hydrosaurus. 
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